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CONVERSION  FACTORS ,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement 

units  as  follows: 

Multiply 

used  in  this  report 

By 

can  be  converted  to  SI  (metric) 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

metres 

gallons  (US  liquid) 

3.785412 

cubic  decimetres 
(litres) 

inches 

25.4 

millimetres 

kips  (force) 

4.448222 

kllonewtons 

kips  (force)  per 

6.894757 

megapascal 8 

square  inch 

megatons  (nuclear 

4.184 

petajoules 

equivalent  of  TNT) 

pounds  (force)  per 

6.894757 

kilopascals 

square  inch 

pounds  (mass) 

0.4535924 

kilograms 

pounds  (mass)  per 

16.01846 

kilograms  per 

cubic  foot  cubic  metre 


NONLINEAR  RESPONSE  OF  CONCRETE  GRAVITY  DAMS  TO 
STRONG  EARTHQUAKE-INDUCED  GROUND  MOTION 

PART  Is  INTRODUCTION 


1.  The  Corps  of  Engineers  (CE)  has  recently  published  (Department  of 
the  Army  (DA)  1985)  criteria  and  guidance  for  assessing  the  seismic  resis- 
tance  of  concrete  gravity  dams.  This  publication  reflects  the  products  of 
research  conducted  by  the  CE  and  others  during  the  last  decade  and  thus 
represents  the  most  up-to-date  information  now  available.  In  contrast  to  the 
static  seismic  coefficient  method  of  analysis  previously  used  (DA  1960)  the 
new  guidance  considers  the  dynamic  properties  of  the  dam,  the  local  seis¬ 
micity  of  the  site  and,  to  some  extent,  the  interaction  among  the  dam, 
reservoir  and  foundation.  However,  the  procedure  reproduced  in  Figure  1  to 
estimate  the  extent  and  significance  of  cracking  experienced  during  extreme 
ground  shaking  is  admittedly  interim  and  lacking  in  extensive  justification. 

2.  The  objective  of  the  research  described  herein  is  to  confirm  or 
appropriately  modify  the  sequence  of  analysis  in  Figure  1.  It  specifically 
encompassed  the  following  contractual  tasks: 

£.  Perform  nonlinear  analyses  on  three  (3)  given  nonoverflow 

gravity  dam  cross  sections,  using  a  sufficient  number  of  strong 
earthquake  ground  motion  records  to  bracket  the  frequency  range 
of  significant  dynamic  response  for  the  gravity  dam  sections. 
The  results  from  these  analyses  will  be  used  to  estimate  the 
extent  of  cracking  possible  and  any  permanent  displacements 
likely  to  occur  along  weak  planes  in  the  foundation,  at  the 
dam- foundation  Interface  or  within  the  dam  cross  section  in 
areas  of  maximum  cracking. 

1).  Evaluate  the  interim  procedure  contained  in  ETL  1110-2-303, 
paragraphs  6  and  7,  pages  1-8  through  1-11,  for  analyzing  the 
response  of  a  gravity  dam  to  a  maximum  credible  earthquake  and 
evaluating  the  results  of  the  analysis,  and  recommend  improve¬ 
ments  based  upon  the  nonlinear  analyses  performed,  as  described 
in  ci  above. 

£.  Develop  a  simplified  method  of  analysis  to  estimate  upper 
bounds  for  permanent  displacements  due  to  vibratory  ground 
motion  in  the  foundation,  at  the  dam-foundation  interface,  and 
within  the  dam  cross  section  without  resorting  to  multiple 
acceleration  time-history  analyses. 


3.  The  remainder  of  this  report  confirms  and  modifies  the  evaluation 
procedure  of  Figure  1.  Sections  2  through  4  describe  our  work  in  accordance 
with  the  foregoing  tasks.  In  Section  5  we  summarize  our  conclusions  and  we 
present  our  recommendations  in  Section  6. 


PART  II:  ADINA  ANALYSES 


4.  The  nonlinear  analyses  of  gravity  dam  response  to  seismic  motion 
were  performed  using  the  general-purpose  finite  element  code  ADINA  (ADINA 
Engineering  1984).  The  cross  section  of  each  structure  was  modeled  using 
9-node ,  Isoparametric,  quadrilateral,  plane  stress  elements  of  unit 
thickness.  The  discretization  of  the  finite  element  meshes  corresponded  to 
SUBROUTINE  INPUTG  of  Cole  and  Cheek  (1985)  wherein  this  degree  of  resolution 
was  found  to  adequately  reproduce  the  stress  distributions  within  gravity 
dams. 

5.  The  constitutive  behavior  of  the  concrete  of  the  dams  was  described 
with  the  ADINA  concrete  material  model  (Bathe  and  Samaswamy  1979).  For  the 
static  and  the  seismic  loadings  imposed,  this  description  provided  an 
essentially  linear  behavior  in  compression  and  a  linear  behavior  in  tension 
up  to  the  stress  level  at  which  cracking  occurred  as  seen  in  Figure  2.  When 
this  level  of  stress  was  reached,  spatial  zones  of  cracked  material  were 
defined.  In  these  zones  the  tensile  stiffness  across  the  cracked  surface  was 
reduced  to  zero  and  the  shear  stiffness  was  reduced  to  half  the  uncracked 
value.  The  effects  of  strain  rate  on  material  behavior  were  approximated  as 
follows.  First,  the  static  stiffness  and  strength  were  related  through 
(American  Concrete  Institute  1977): 
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(1) 


in  which: 

-  Initial  modulus  of  elasticity,  psi. 
f^  “  Uniaxial  compressive  strength,  psi. 
w  -  Unit  weight,  150  pcf. 

Next,  the  static  uniaxial  tensile  strength  in  psi,  f^  ,  was  found  from 
(Raphael  1982): 

f  -  1.7  f,2/3  (2) 

t  c 


6 


Static  uniaxial  ultimate  compressive  stresses  f  were  estimated  from 
typical  stress-strain  curves  corresponding  to  f^  (Winter  and  Nilson 
1979).  In  accordance  with  recent  dynamic  biaxial  test  (Mlakar,  Vltaya-Udom, 
and  Cole  1984),  a  dynamic  strength  amplification  was  estimated  from: 
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(3) 


in  which  one  fourth  of  the  dam's  fundamental  period  was  substituted  for  the 
rise  time  of  loading  in  msec,  tf  .  The  constant  strength  parameters  used 
in  the  ADINA  calculations  were  then: 


a  -  a  f  (4) 

a  =  a  f  (5) 

c  c 


a  =  a  f  (6) 

u  u 


In  which: 

*  Uniaxial  cut  off  tensile  stress. 
ac  =*  Uniaxial  maximum  compressive  stress. 

8u  =  Uniaxial  ultimate  compressive  stress. 

The  uniaxial  compressive  strain  corresponding  to  3^  and  the  uniaxial 

ultimate  compressive  strain  e^  were  estimated  from  the  static  stress-strain 
curves  corresponding  to  (Winter  and  Nilson  1979).  This  is  supported 

by  the  finding  of  Mlakar,  Vltaya-Udom,  and  Cole  (1984)  that  the  strain  at 
concrete  failure  is  independent  of  rate  effects.  Finally,  the  tangent 
modulus  at  zero  strain  Eq  was  approximated  to  be 
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(7) 


for  consistency  with  the  foregoing.  The  numerical  values  of  these  material 
parameters  Input  to  ADINA  are  listed  for  each  dam  in  Table  1.  In  each  case, 
a  Poisson's  ratio  of  0.2  was  used. 


6.  In  Che  nonlinear  analysis,  both  static  and  seismic  loadings  were 
considered.  The  self-weight  of  the  cross  section  was  applied  as  a  mass 
proportional  body  force.  The  hydrostatic  loading  of  the  reservoir  was 
accomplished  through  a  set  of  pressure  loadings  along  the  upstream  surface  of 
the  structure.  Both  of  these  loadings  were  linearly  increased  from  zero 
initially  to  their  full  static  value  in  a  rise  time  of  2  sec.  This  rise  time 
was  selected  to  achieve  equilibrium  under  the  static  loading  prior  to  the 
arrival  of  the  strong  earthquake  motion.  The  horizontal  and  the  vertical 
earthquake  loading  were  applied  as  uniformly  prescribed  displacements  along 
the  base  of  the  dams.  Finally,  the  hydrodynamic  loading  of  the  reservoir  was 
approximated  by  adding  concentrated  nodal  masses  on  the  upstream  face 
corresponding  to  the  distribution  (Chopra  1978) 


m  (z) 
a 


p.  (z  ,  m  ) 
1 _ s 

v(z) 


(8) 


in  which: 


z  -  Vertical  coordinate. 

gp^  -  Hydrodynamic  pressure  in  fundamental  mode. 

uj  ■  Fundamental  frequency  including  hydrodynamic  effects. 
3 

■  Fundamental  mode  shape, 
g  -  Acceleration  of  gravity. 


7.  The  foregoing  constitutes  a  problem  which  is  both  materially 
nonlinear  and  dynamic.  It  was  solved  through  implicit  integration  with 
respect  to  time  using  the  trapezoidal  rule  with  a  step  of  0.02  sec.  This 
value  was  the  sampling  interval  of  the  recorded  earthquake  loadings  and  is 
thus  presumably  smalL  enough  to  describe  the  energy  of  the  ground  motion. 
Consequently,  the  response  of  structures  to  this  excitation  is  also 
adequately  represented  with  this  step.  The  BFGS  matrix  update  method  (Bathe 
1982)  was  employed  to  iteratively  solve  the  equilibrium  equations  with  an 
energy  convergence  criteria.  Stiffness  reformation  and  equilibrium  iteration 
were  performed  in  each  time  step. 

8.  Three  dams  were  analyzed  in  this  study.  They  were  selected  to 
characterize  that  segment  of  the  Corps'  population  of  such  structures  which 
is  of  the  greatest  seismic  interest.  Two  actual  projects  and  a  standard 
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cross  section  representative  of  other  structures  are  included.  The  three 
structures  considered  are  a  subset  of  those  linearly  analyzed  in  Cole  and 
Cheek  (1985)  so  as  to  facilitate  direct  comparisons  with  this  work.  The 
geometric  parameters  describing  these  structures  are  compared  in  Table  2. 

9.  The  shortest  structure  is  the  Richard  B.  Russell  Dam  of  the 
Savannah  District.  This  project  is  185  ft  tall  and  has  a  modulus  of 
elasticity  equal  to  3  million  psi.  This  dam  was  included  because  a  number  of 
previous  studies  have  been  conducted  on  it  (Norman  1979)  (Norman  and  Stone 
1979)  (Chiarito  and  Mlakar  1983).  The  finite  element  grid  used  in  the 
analysis  is  shown  in  Figure  3.  It  contains  65  elements  and  297  nodes  having 
594  degrees  of  freedom  in  all. 

10.  The  second  structure  is  a  "standard"  dam  having  a  horizontal 
foundation,  a  near  vertical  upstream  face,  a  rectangular  crown  and  a 
downstream  face  with  a  slope  of  10  on  12  as  seen  in  Figure  4.  The  structural 
height  is  300  ft  which  is  approximately  the  ninetieth  percentile  of  the 
Corps'  population.  Computations  were  made  using  a  modulus  of  concrete  of 

3  million  psi  which  is  representative  of  Corps'  construction  practice.  The 
finite  element  mesh  of  this  structure  has  594  degress  of  freedom  in  all. 

11.  The  tallest  dam  considered  is  the  Dworshack  Dam  on  the  Clearwater 
River  in  Idaho.  This  structure  is  638  ft  high,  and  is  made  of  concrete  with 
an  estimated  5  million  psi  modulus.  It  is  included  to  be  representative  of 
the  very  tallest  concrete  dams  owned  by  the  CE.  The  finite  element 
discretization  of  this  dam  has  80  elements,  363  nodes,  726  degress  of  freedom 
and  is  shown  in  Figure  5. 

12.  Preliminarily  to  the  nonlinear  analyses,  the  first  four  periods 
and  mode  shapes  of  each  structure  without  reservoir  effects  were  found.  The 
results  in  Table  3  and  Figures  6,  7  and  8  are  consistent  with  previous 
studies  of  gravity  dam  dynamics  (Chopra  and  Chakrabarti  1971)  (Norman  1979) 
(Norman  and  Stone  1979). 

13.  The  earthquake  loading,  Figure  9,  used  in  the  analyses  was  the 
N65W  and  vertical  components  of  the  1966  Parkfield,  California  earthquake 
recorded  at  Temblor  No.  2  Station,  CIT  File  Nos.  B037-1  and  B037-3.  This 
loading  is  representative  of  strong  earthquake  motions  likely  to  be  encoun¬ 
tered  at  rock  sites  upon  which  all  CE  concrete  dams  are  founded.  Table  4 
lists  three  pertinent  parameters  of  22  earthquake  recordings  at  such  sites  as 


tabulated  by  Lai  (1980).  Mote  that  the  peak  ground  acceleration  A  of 

III 

the  selected  loading  Is  one  of  the  greatest  in  Table  4 ,  and  thus  most  likely 
to  significantly  excite  the  nonlinear  response  of  concrete  dams.  The 

parameter  o>  characterizes  the  ground  frequency  of  the  power  spectral 

© 

density  function  (Kanai  1961)  (Tajimi  1960).  Mote  in  Figure  10  that  this 
parameter  is  independent  of  the  peak  ground  acceleration  and  that  the  value 
for  the  selected  earthquake  is  representative  of  the  entire  set  of  records. 
The  equivalent  strong  motion  duration  Sq  ,  as  defined  by  Vanmarcke  and  Lai 
(1980)  is  also  listed  in  Table  4.  In  Figure  11,  note  that  this  parameter 
tends  to  decrease  with  increasing  peak  ground  acceleration.  The  record 
selected  is  seen  to  be  representative  of  the  stronger  recordings  in  this 
Important  respect.  It  should  also  be  noted  that  the  Parkfleld  loading  chosen 
for  this  study  is  one  of  those  used  by  Cole  and  Cheek  (1985),  and  thereby 
facilitates  direct  comparison  with  their  work. 

14.  In  this  section  the  structural  geometry,  material  behavior, 
loading  conditions  and  nonlinear  solution  scheme  of  the  ADINA  analyses  have 
been  described.  The  results  of  these  calculations  are  presented  in  the 
following  section.  Therein,  these  results  are  also  used  to  judge  the 
structural  evaluation  procedure  advocated  in  DA  (1985). 
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PART  III:  EVALUATION  OF  ETL  1110-2-303 


15.  In  the  previous  section,  we  described  the  manner  In  which  the 
ADINA  finite  element  code  was  used  to  model  the  nonlinear  response  of 
concrete  gravity  dams  to  strong  earthquake- Induced  ground  motion.  Herein, 
the  results  of  these  analyses  are  discussed  and  used  to  evaluate  the  Corps' 
criteria  and  guidance  for  assessing  the  seismic  resistance  of  such  structures 
(DA  1985).  The  remainder  of  this  section  is  organized  by  each  of  the  four 
dam-earthquake  combinations  which  were  analyzed. 

Russell  Dam 


16.  The  first  dam  to  be  discussed  is  the  Russell  Dam.  When  It  was 
subjected  to  the  Parkfleld  motion  with  the  ADINA  model,  no  cracked  zones 
propagated  through  any  portion  of  the  cross  section.  This  case  was  thus  of 
trivial  nonlinear  Interest.  It  does,  however,  provide  an  Interesting 
instance  for  the  evaluation  of  ETL  1110-2-303. 

17 .  The  results  of  all  four  dam-earthquake  analyses  with  the  guidance 

of  this  ETL  are  summarized  In  Table  5.  Therein,  05  *  07  ’  and  fo  10 

are  the  maximum  stresses  calculated  with  the  simplified  procedure  (Cole  and 
Cheek  1985)  for  5,  7,  and  10  percent  of  critical  damping,  respectively.  The 
guidance.  Figure  1,  first  prescribes  an  analysis  with  5  percent  damping.  The 
combined  maximum  stress  for  Russell  Dam  is  640  psi  at  the  upstream  edge  of 
the  base.  As  this  exceeds  0.20  f^  ,  a  reanalysis  with  10  percent  damping 

is  performed.  The  maximum  stress  is  reduced  to  530  psi  at  the  upstream  edge 
of  the  base.  Since  this  value  exceeds  0.10  f^  ,  a  sliding  stability 
analysis  on  the  base  plane  is  required. 

18.  In  this  case,  the  ETL  guidance  is  conservative  when  judged  by  the 
ADINA  results.  The  ETL  requires  a  sliding  stability  analysis  along  a  hypo¬ 
thetical  cracked  plane  while  the  ADINA  analyses  indicates  no  propagation  of 
cracked  zones  through  the  dam.  There  are  two  reasons  for  this  difference. 
First,  the  response  of  this  dam  to  the  equivalent  static  loading  of  the  ETL 
overestimates  the  dynamic  response  for  this  earthquake  (Cole  and  Cheek 
1985).  Secondly,  the  tensile  strength  modeled  in  ADINA,  8  -  380  psi  in 
Table  1,  accounts  for  some  dynamic  Increase  above  the  0.15  f^  of  the  ETL. 
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19.  To  obtain  a  significant  nonlinear  response  in  this  structure,  it 
was  also  subjected  to  a  base  motion  of  three  times  the  amplitude  of  the 
recorded  Parkfleld  Earthquake.  In  the  ADINA  analysis  of  this  case,  a  cracked 
zone  does  propagate  through  the  cross  section.  An  indication  of  this 
propagation  is  obtained  from  Figure  12  which  shows  the  total  number  of 
elemental  integration  points  at  which  cracking  has  occurred  as  a  function  of 
time.  Figure  13  shows  the  extent  of  the  cracked  zones  at  the  key  instants  of 
time  identified  in  Figure  12.  No  cracking  occurs  prior  to  the  arrival  of  the 
strong  ground  shaking  as  one  would  expect.  Cracking  initiates  near  the 
upstream  edge  of  the  base  and  temporarily  stabilizes  until  3.84  sec. 
Thereafter,  it  propagates  toward  the  downstream  edge.  Simultaneously,  a 
cracked  zone  initiates  near  the  change  in  upstream  slope.  This  zone 
stabilizes  while  the  one  at  the  base  propagates  through  the  downstream  face. 

20.  The  analysis  of  this  situation  prescribed  by  the  ETL  is  reported 
in  the  third  column  of  Table  5.  In  this  case,  the  ETL  guidance  leads  to 
results  which  are  consistent  with  the  ADINA  analysis  in  the  essential 
features.  The  maximum  stresses  computed  by  the  simplified  method  occur  at 
the  upstream  edge  of  the  base  and  the  criteria  require  a  sliding  stability 
analysis  along  this  horizontal  plane.  The  more  comprehensive  ADINA  analysis 
indicates  cracking  from  the  upstream  face  to  the  downstream  face  at  this  same 
elevation. 


Standard  Dam 


21.  The  second  structure  to  be  considered  is  the  300-ft-tall  Standard 
Dam.  In  the  ADINA  analysis  with  the  recorded  Parkfleld  motion,  two  cracked 
zones  propagate  through  the  cross  section.  The  progression  of  cracking  with 
time  is  shown  in  Figure  14  and  the  zones  of  cracking  at  critical  times  are 
presented  in  Figure  15.  Note  that  a  cracked  zone  initiates  at  the  upstream 
edge  of  the  base,  propagates  toward  the  downstream  edge,  but  stabilizes 
before  transecting  the  cross  section.  Subsequently,  two  cracked  zones  form 
on  the  downstream  surface  at  the  elevation  of  slope  change  and  at  a  slightly 
lower  elevation.  These  then  quickly  propagate  to  the  upstream  surface. 
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22.  The  ETL  analysis  of  the  corresponding  case  Is  summarized  In  the 
fourth  column  of  Table  S.  In  this  case,  the  ETL  guidance  Is  In  partial 
agreement  with  the  ADINA  model.  The  ETL  procedure  does  Indicate  that  this 
structure  will  crack  through  the  cross  section.  However,  the  ETL  procedure 
predicts  this  cracking  to  occur  at  the  base  while  the  ADINA  model  shows  it  to 
occur  in  the  upper  portion  of  the  dam. 

Dworshak  Dam 


23.  The  final  dam  considered  is  the  638-ft-tall  Dworshak  structure. 

In  the  ADINA  analysis,  a  cracked  zone  transects  the  cross  section  at  the 
elevation  of  change  in  downstream  slope.  The  progression  of  cracking  in  time 
and  space  is  shown  in  Figures  16  and  17.  Cracking  initiates  at  the  upstream 
edge  of  the  base  and  soon  stabilizes.  Later,  a  cracked  zone  propagates  from 
the  downstream  face  at  the  elevation  of  slope  change  and  another  propagates 
from  a  lower  elevation  on  this  face.  The  latter  propagates  toward  a  fourth 
zone  which  subsequently  forms  on  the  upstream  face.  However,  the  transection 
of  the  cross  section  occurs  in  the  zone  at  the  elevation  of  downstream  slope 
change . 

24.  The  ETL  analysis  of  Dworshak  Dam  is  summarized  in  the  last  column 
of  Table  5.  These  results  are  partially  consistent  with  those  of  the  ADINA 
calculation.  The  ETL  procedure  indicates  that  cracking  is  expected  through 
the  base  of  the  structure.  The  ADINA  model  also  indicates  cracking  through 
the  structure  but  near  the  elevation  of  change  in  downstream  slope. 

Discussion 


25.  When  excited  by  an  earthquake  motion  strong  enough  to  initiate 
cracking,  the  dams  considered  cracked  completely  through  their  cross 
sections.  This  is  consistent  with  limited  scale  model  tests  (Niwa  and  Clough 
1980).  It  is  noted  here  that  in  some  cases  the  cracked  zones  grew  through 
stable  stages.  In  other  cases,  the  cracking  was  virtually  instantaneous  once 
it  began.  The  practical  implication  of  this  may  be  that  dams  which  crack 
during  an  earthquake  do  crack  completely  through  their  cross  section. 


However,  further  analytical  and  experimental  study  Is  required  to  support 
this  conclusion  in  general. 

26.  The  three  dams  considered  cracked  completely  at  various  eleva¬ 
tions.  It  is  noted  that  the  shortest  structure  cracked  through  its  base 
while  the  taller  ones  cracked  through  the  upper  elevations.  The  elevation  of 
complete  cracking  may  be  sensitive  to  the  flexibility  of  the  foundation  which 
was  not  considered  in  this  nonlinear  analysis.  However,  a  state-of-the-art 
linear  analysis  which  included  this  feature  (Fenves  and  Chopra  1984) 
indicates  that  the  maximum  total  stress  occurs  at  the  upstream  edge  of  the 
base.  This  is  the  samr  location  at  which  cracking  Initiated  in  each  of  the 
nonlinear  analyses  conducted  herein. 

27.  The  simple  procedure  of  ETL  1110-2-303  conservatively  evaluated 
the  safety  of  each  of  the  four  dam-earthquake  combinations  considered 
comprehensively  in  this  work.  But  in  two  cases,  the  simple  linear  procedure 
was  misleading  in  that  the  elevation  of  cracking  through  the  cross  section 
differed  from  that  observed  in  the  nonlinear  analysis.  It  is  Important  that 
this  elevation  be  correctly  predicted  so  that  a  subsequent  estimate  of  perma¬ 
nent  displacement  along  the  cracked  plane  is  meaningful.  This  prediction 
requires  that  the  inelastic  behavior  following  initial  cracking  be  incorpo¬ 
rated  in  the  evaluation  procedure.  This  might  be  simply  achieved  with  a 
static  nonlinear  resistance  function  and  an  inelastic  response  spectrum. 
Further  study  is  suggested  to  revise  the  procedure  in  this  way. 

Summary 


28.  In  this  section  we  have  discussed  the  results  of  ADINA  analyses  of 
four  dam-earthquake  combinations.  The  location  and  extent  of  cracking 
calculated  by  ADINA  have  further  been  used  to  evaluate  the  procedures 
recommended  in  ETL  1110-2-303.  In  the  following  section,  a  procedure  Is 
developed  to  estimate  the  permanent  relative  displacement  likely  to  occur 
along  cracks  in  these  zones. 


14 


PART  IV:  PERMANENT  DISPLACEMENTS 


29.  In  Che  foregoing  section  It  was  seen  that  extremely  strong  ground 
shaking  can  cause  cracking  through  the  cross  section  of  a  concrete  gravity 
dam.  In  fact,  this  may  have  occurred  at  the  Royna  Dam  during  the  December 
1967  earthquake  (Chopra  and  Chakrabarti  1971).  The  Corps'  analytical 
procedure  for  the  maximum  credible  earthquake  (Figure  1),  requires  an 
estimate  of  permanent  displacement  in  such  a  case.  A  simple  procedure  for 
this  estimation  is  developed  and  Illustrated  in  this  section. 

30.  In  the  proposed  procedure,  the  cracked  dam  is  idealised,  as  shown 

in  Figure  18.  The  cracked  surface  is  assumed  to  be  planar  and  inclined  at 
angle  6  to  the  hr  'zontal.  The  portion  of  the  monolith  above  this  crack  is 
assumed  to  be  a  ~  N)dy  subject  to  sliding  along  the  crack  plane  as  the 
portion  of  oelow  the  crack  moves  horizontally  due  to  the  ground 

motir  at ion  resembles  the  sliding  block  analysis  widely 

emplo,  .a  embankments  (Newmark  1965)  (Srrma  1975)  (Franklin  and 

Chang  ,  (Hynes-Griffin  and  Franklin  1984)  and  should  be  no  less 
applicable  to  the  problem  at  hand. 

31.  The  general  application  of  the  sliding  block  analyses  has  been 

facilitated  by  a  recent  probabilistic  treatment  (Lin  and  Whitman  1986).  In 

Figure  19,  the  expected  permanent  displacement  D  is  presented  as  a  function 

of  peak  base  acceleration  A  and  the  limiting  acceleration  A  •  The 

max  c 

result  for  rock  sites  is  applicable  as  all  CE  gravity  dams  are  so  founded. 
Observe  that  no  permanent  displacement  is  expected  if  A£  >  A  .  For 
Ac  ^  Amax  '  the  ProbaWllty*  p  »  of  exceeding  various  displacement 
thresholds,  Dq  ,  can  subsequently  be  obtained  from  Figure  20. 

32.  For  the  cracked  gravity  dam,  the  limiting  acceleration  A£  at 
which  sliding  Impends  follows  from  equilibrium  of  Figure  18,  wherein: 

W  -  Weight  of  dam  above  crack. 

H  -  Hydrostatic  force  on  dam  above  crack. 

N  -  Normal  force  across  crack. 

F  ■  Frictional  force  along  crack. 

A  "  Base  horizontal  acceleration. 


For  a  coefficient  of  friction  p  end  e  gravitational  acceleration  g  ,  the 
Uniting  value  for  downstream  eliding  ia: 


m  yi  cos  8  -  sin  6  H 

g  p  sin  8  +  coa  6  W 


(9) 


The  Uniting  value  correaponding  to  upatrean  eliding  ia  given  by: 

^  ,  tLcos  9  t  sin  6  +  H 
g  cos  6  -  psin  6  W 

and  thla  nay  be  aore  critical  for  anall  hydrostatic  forces  and  cracks  whose 
normal  la  inclined  upatrean  (negative  values  of  8).  This  is  seen  in  the 
plot  of  Uniting  acceleration  aa  a  function  of  crack  inclination  in  Figure  21 
for  u  -  1.0  and  H/W  *0.1  .  Observe  that  A£  decreases  for  sharply 
inclined  cracks.  As  H/W  increases,  the  upstream  critical  acceleration 
increases  while  the  downstrean  value  decreases  from  that  shown.  For  increas¬ 
ing  p  ,  both  upstream  and  downstream  A  Increase. 

c 

33.  Sone  indication  of  an  appropriate  value  of  the  frictional 
coefficient  p  la  obtained  from  direct  shear  tests  of  concrete  and  rock 
speclnens  which  are  sunmarlzed  in  Table  6.  The  values  for  the  precut 
specimens  undoubtedly  underestimate  the  resistance  provided  by  naturally 
cracked  concrete.  However,  the  consistency  between  precut  concrete  and  rock 
frictional  coefficients  suggests  that  the  coefficients  for  natural  rock 
joints  nay  also  represent  those  for  natural  concrete  cracks.  From  model 
tests  on  concrete  gravity  dams  (Niwa  and  Clough  1980),  we  anticipate  any 
earthquake- Induced  cracked  surfaces  to  be  relatively  rough.  From  the  limited 
information  available,  a  value  of  p  ■  1.0  is  thus  preliminarily  proposed 
for  use  in  the  sliding  block  analysis.  However,  a  further  investigation  of 
this  Important  parameter  nay  be  desirable. 

34.  The  peak  base  acceleration  A  should  be  chosen  to  reflect 

max 

the  amplification  of  ground  motion  at  the  level  of  cracking  in  the  structure. 
The  absolute  acceleration  a  at  elevation  z  and  time  t  is  (Biggs  1964): 

a(z,  t)  -  a  (t)  +  ?  r  il°(t)  (z)  (11) 

g  n  «•  1  n  n  n 
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wherein: 


a  (c)  ■  Ground  acceleration  function, 
g 

n  ■  Mode  number. 

r  ”  Modal  participation  factor, 
n 

<P  (z)  -  Mode  shape  function, 
n 

ii°(t)  ■  Relative  acceleration  response  of  slngle-degree-of-freedom 
system  having  period  and  damping  of  mode  n  . 


If  only  the  first  mode  of  vibration  is  considered,  a  reasonable  root  mean 
square  estimate  of  peak  base  acceleration  at  the  elevation  of  cracking  is 
(Newmark  and  Rosenbleuth  1971): 


Amax  ‘  [>  +  (fl  Sal  ^  <12> 


in  which  a _  Is  the  maximum  ground  acceleration  and  S  ,  Is  the 

gm  al 

ordinate  of  the  psuedoacceleratlon  response  spectrum  at  the  fundamental 
period  corresponding  to  a  (t)  .  Consistency  with  the  simplified  stress 

O 

analysis  of  gravity  dams  (Chopra  1978)  Is  achieved  by  taking  -  4  and 


2 

^(z)  -  -0.0002  +  0.1427  (f-)  +  0.6683  (jp) 

s  s 

3  4 

-0.8508  (^f-)  +  1.0399  (~) 


(13) 


where  Is  the  structural  height  of  the  dam.  The  amplification  of 

ground  motion  Implied  by  Equation  12  Is  shown  in  Figure  22.  Note  that 

A  Increases  with  the  elevation  of  cracking  and  with  the  spectral 
max 

acceleration  as  one  would  expect. 

35.  To  illustrate  this  procedure,  let  us  examine  the  experience  of  the 
Koyna  Dam  during  the  strong  1967  earthquake.  From  the  evidence  of  upstream 
and  downstream  surflcial  cracking  (Chopra  and  Chakrabarti  1971),  we  assume  a 
horizontal  crack,  as  shown  in  Figure  23,  where  elevations  are  referenced  in 
feet  to  the  Koyna  Reduced  Level  (KRL).  For  a  concrete  density  of  165  lb/ft^ 


(Saini  and  Krishna  1974),  H/W  -  0.19  .  With  y  -  1.0  ,  Equation  9  then 
indicates  the  Uniting  acceleration  to  be  A£  ■  0.81  g  .  For  the  geonetry 
of  Figure  23  and  E  -  3.6  x  10^  psi  (Saini  and  Krishna  1974),  the  funda¬ 
mental  period  Including  the  influence  of  the  stored  water  is  T  *  0.30  sec 
(Chopra  1978).  Corresponding  to  the  transverse  motion  recorded  on 

11  December  1967  (Chopra  and  Chalcrabarti  1971),  we  have  a  *  0.49  g  and 

g® 

S  ■  0.56  g  for  5  percent  of  critical  damping.  Equation  12  then 

estimates  a  peak  base  motion  A  ■  0.87  g  .  For  A  /A  »  0.93  , 

max  c  max 

Figure  19  implies  an  expected  permanent  displacement  D  <  0.1  x0.87  = 

0.087  cm.  As  no  offset  in  the  cracked  monoliths  was  observed  following  the 
1967  earthquake,  we  gain  confidence  in  the  procedure  developed  herein.  It  is 
also  noted  that  ^max  “  0.87  g  estimated  with  Equation  12  compares 
favorably  with  0.8  g  calculated  through  numerical  integration  of  the 
equations  governing  transverse  vibration  (Saini  and  Krishna  1974). 

36.  As  a  further  application  of  the  method,  consider  the  Russell  Dam 

subjected  to  the  Parkfield  Earthquake  motion.  As  shown  in  Figure  24,  we 

hypothetically  assume  that  a  horizontal  crack  occurs  at  the  elevation  of 

change  in  downstream  slope  where  the  stress  is  a  local  maximum  in  the 

simplified  linear  analysis  (Cole  and  Cheek  1985).  For  the  conditions  shown, 

the  limiting  acceleration  Ac  -  0.95  g  follows  from  Equation  8. 

Equation  12  indicates  a  peak  acceleration  A  -  1.77  g  .  The  expected 

max 

permanent  displacement  D  -  4.4  cm  then  follows  from  Figure  19.  Figure  20 
may  be  used  to  quantify  the  uncertainty  in  this  estimate  due  to  the  random¬ 
ness  of  ground  motion.  For  example,  an  exceedance  probability  of  P  -  0.1 

corresponds  to  the  permanent  displacement  D  •  2.2  x  4.4  «  9.8  cm  .  This 

o 

is  a  conditional  probability  given  the  values  of  limiting  and  peak 

accelerations.  As  considerable  uncertainties  exist  with  respect  to  A^ 

and  A  ,  the  permanent  displacement  corresponding  to  an  unconditional 
max 

probability  of  exceedance  may  differ  from  that  given  by  Figure  20.  It  should 
be  noted  that  these  displacements  are  conservative  as  the  nonlinear  analyses 
of  the  previous  section  did  not  Indicate  any  cracking  in  the  upper  portion  of 
this  dam  for  the  Parkfield  motion. 

37.  The  sliding  block  analysis  can  thus  be  applied  to  simply  estimate 
the  permanent  displacements  of  gravity  dams  cracked  by  strong  ground  motion. 
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In  the  analysis,  the  Uniting  acceleration  is  calculated  from  equilibrium  of 
the  monolith  above  the  crack  and  the  amplification  of  ground  motion  is 
computed  from  an  approximation  consistent  with  the  simplified  stress  analysis 
of  gravity  dams.  The  procedure  leads  to  a  conservative  result  in  each  of  the 
examples  Illustrated.  The  method  may  be  improved  through  an  experimental 
examination  of  the  frictional  coefficient  employed  to  compute  the  limiting 
acceleration.  An  analytical  study  of  the  amplitude  and  frequency  modulation 
of  ground  motion  within  the  cracked  structure  may  also  lead  to  a  better 
procedure . 


PART  V:  CONCLUDING  REMARKS 


38.  The  funding  for  this  contract  permitted  the  investigation  of  a 
rather  limited  number  of  cases.  Thus,  broadly  applicable  conclusions  are  not 
supported  by  this  study  alone.  However,  the  following  observations  can  be 
made: 

a.  When  excited  by  an  earthquake  motion  of  sufficient  strength  to 
initiate  cracking,  the  three  dams  considered  cracked  completely 
through  their  cross  sections.  In  some  cases,  this  cracking 
progressed  through  stable  stages  and  in  others,  it  was  virtually 
instantaneous. 

b.  The  cracking  through  the  dam  cross  sections  occurred  at  the  base, 
at  the  elevation  of  downstream  slope  change,  and  at  a  lower 
elevation.  The  shortest  structure  considered  cracked  through  its 
base  while  the  taller  ones  cracked  through  the  upper  elevations. 

£.  The  recently  published  procedure  of  ETL  1110-2-303  conservatively 
evaluated  the  safety  of  the  four  cases  considered.  However,  in 
two  cases  the  simplified  analysis  procedure  incorrectly  located 
the  elevation  at  which  cracking  occurs  through  the  cross  section. 

<1.  The  sliding  block  analysis  is  a  rational  basis  to  estimate  the 
permanent  relative  displacements  along  cracked  planes  through  a 
dam  cross  section. 


PART  Vis  RECOMMENDATIONS  FOR  FURTHER  STUDY 


39,  The  scope  of  this  Investigation  was  Halted  to  first  learn  what 
might  be  achieved  before  undertaking  a  more  comprehensive  examination  of  the 
nonlinear  response  of  concrete  gravity  dams  to  strong  earthquake-induced 
ground  motion.  Based  on  the  observations  of  the  preceding  section,  the 
following  specific  recommendations  for  further  Investigation  are  offered: 

£.  Further  dynamic,  nonlinear,  finite  element  analyses  should  be 
conducted  of  earthquake- loaded  gravity  dams  to  genrallze  the 
tentative  observations  about  the  location  and  the  extent  of 
cracking.  In  these  analyses,  additional  earthquake  motions 
recorded  on  rock  should  be  used  to  excite  the  structure.  Also,  a 
simple  model  of  the  most  essential  effects  of  foundation 
flexibility  should  be  incorporated. 

t>.  The  simplified  method  of  analysis  should  be  revised  so  that  It 
correctly  locates  the  elevation  at  which  cracking  occurs  through 
the  cross  section  of  gravity  dams.  This  can  probably  be  achieved 
by  including  the  most  essential  features  of  foundation  compliance 
and  nonlinear  postcracking  behavior. 

£.  The  frictional  resistance  along  naturally  cracked  concrete 

surfaces  should  be  experimentally  studied  to  obtain  appropriate 
frictional  coefficients  for  the  sliding  block  analysis. 
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Table  1 

Parameter  Values  Used  In  ADINA  Concrete  Model 


Parameter 

Russell  Dam 

Standard  Dam 

Dworshak  Dam 

Eq  ,  106  psl 

3.34 

3.33 

5.52 

3t  .  P«i 

380 

380 

740 

3c  ,  psi 

-3040 

-3010 

-8300 

e 

c 

-0.002 

-0.002 

-0.002 

-2230 

-2220 

-3050 

e 

-0.004 

-0.004 

-0.003 

u 


Table  2 

Geometric  Parameter  Values  for  Dams  Analyzed 


Parameter 

Russell  Dam 

Standard  Dam 

Dworshak  Dam 

Dam  height,  ft 

185 

300 

638 

Upstream  slope 

1/12 

1/12 

0 

Upstream  slope  height,  ft 

143 

250 

— 

Crest  width,  ft 

17 

21 

27 

Downstream  slope 

2/3 

7/10 

4/5 

Downstream  slope  height,  ft 

160 

270 

585 

Base  width,  ft 

136 

231 

495 

Reservoir  height,  ft 

170 

285 

570 

Table  3 

Natural  Periods  la  Seconds  of  Dams  Analyzed 


Mode  Russell  Dam  Standard  Dam  Dworshak  Dam 

1  0.1581  0.2375  0.3676 

2  0.0743  0.1129  0.1849 


Table  4 

Parameters  of  Earthquakes  Recorded  at  Rock  Sites  (Lai  1980) 


err  no. 

A« 

“g 

rad/ sec 

S 

o 

sec 

A015-1 

0.083 

33.50 

1.34 

A015-2 

0.105 

37.45 

2.02 

B037-1* 

0.269 

24.61 

1.39 

B037-2 

0.347 

21.21 

0.65 

B038-1 

1 

0.014 

25.11 

8.27 

B038-2 

0.012 

29.73 

8.07 

B040-1 

0.041 

20.47 

7.01 

B040-2 

0.046 

21.22 

5.65 

D056-1 

0.315 

23.96 

2.71 

D056-2 

0.271 

17.29 

7.07 

V314-1 

0.064 

7.88 

21.62 

V314-2 

0.097 

6.54 

13.48 

V319-1 

0.054 

21.49 

3.40 

V319-2 

0.036 

17.11 

7.05 

V331-1 

0.041 

33.06 

1.63 

V331-2 

0.037 

45.80 

1.27 

W334-1 

0.142 

18.12 

2.42 

W334-2 

0.198 

25.02 

1.26 

W335-1 

0.071 

42.21 

2.98 

W335-2 

0.056 

45.99 

2.26 

W336-1 

0.057 

35.99 

3.36 

W336-2 

0.071 

32.96 

1.90 

Mean 

0.110 

26.67 

4.86 

4 


*  Horizontal  notion  used  In  analyses. 
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Table  6 

Direct  Shear  Testa  of  Precut  Concrete  and  Rock  Specimens 


Condition 


Residual 


Reference 


Concrete  on  dolomite,  precut 
Dolomite  on  dolomite,  precut 
Natural,  smooth  dolomite  joint 

Concrete  on  dolomite,  precut 
Dolomite  on  dolomite,  precut 
Natural,  medium  rough  dolomite  joint 

Concrete  on  very  hard  sandstone,  precut 
Very  hard  sandstone  on  sandstone,  precut 
Natural  very  hard  sandstone  joint 
Hard  sandstone  on  sandstone,  precut 
Natural  hard  sandstone  joint 
Shaley  sandstone  on  sandstone,  precut 


0.49 

Stowe  1978 

0.40 

Stowe  1978 

0.42 

Stowe  1978 

0.49 

Stowe  et 

al.  1980 

0.40 

Stowe  et 

al.  1980 

1.04 

Stowe  et 

al.  1980 

0.63 

Thornton 

et  al.  1981 

0.65 

Thornton 

et  al.  1981 

0.62 

Thornton 

et  al.  1981 

0.50 

Thornton 

et  al.  1981 

1.17 

Thornton 

et  al.  1981 

0.62 

Thornton 

et  al.  1981 

Figure  la 


Figure  1 .  Sequence  of  analysis  for  the  maximum  credible 
earthquake  (DA  1985)  (Sheet  1  of  2) 
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Figure  5.  Finite  element  grid  of  Dworshak  Dam 
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Figure  6.  Characteristic  shapes  of  Russell  Dam 
(Sheet  1  of  4) 
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Figure  7.  (Sheet  2  of  4) 


500.00.0 


500.0  1000.0  1500.0  2000.0  2500.0  3000.0 

Y 


d.  Mode  4 


figure  7.  (Sheet  4  of  4) 
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Figure  8.  (Sheet  2  of  4) 
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Figure  12.  Cracking  of  Russell  Dam  with  tripled  Parkfield  lo  iding 


a.  t  -  3.84  sec 


Figure  13.  Cracked  zones  of  Russell  Dam  with  tripled 
Parkfield  loading  (Sheet  1  of  4) 


b.  t  *  3.92  sec 


Figure  14.  Cracking  of  Standard  Dam 


a.  t  ■  4.00  sec 


Figure  15.  Cracked  zones  of  Standard  Dam  (Sheet  1  of  4) 


b.  t  =  4.40  sec 


Figure  16.  Cracking  of  Dworshak  Dam 


a.  t  *  4.04  sec 

Figure  17.  Cracked  zones  of  Dworshak  Dam 
(Sheet  1  of  4) 


b .  t  *  4.16  sec 
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Figure  17.  (Sheet  3  of  4) 
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Figure  18.  Free  body  dlagraa  of  cracked  daa 
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Figure  19.  Normalized  ospscted  permanent  displacements  for  sliding 
blocks  on  different  soil  sites  (Lin  and  Whitman  1986) 
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Figure  20.  Normalized  conditional  distributions  of  permanent  dis¬ 
placements  for  sliding  blocks  (Lin  and  Whitman  1986) 


Figure  21.  Critical  acceleration  of  cracked  dam  for  y  -  1.0  and  H/W 
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